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SUMMARY 

Some  aspects  of  the  theory  of  a  gas  laser  recently 
developed  by  W.  E.  Lamb,  Jr.  are  recast  in  a  form  which  more 
fully  displays  the  role  played  by  the  particle  dynamics.  The 
Wigner  distribution  function  is  used  to  derive  kinetic 
equations  which  govern  the  external  center  of  mass  motion  of 
the  two-level  systems  as  well  as  their  internal  dynamics. 

The  effect  of  long  range  forces  is  discussed  by  treating  the 
collision  integral,  in  a  manner  similar  to  that  employed  in 
plasma  kinetic  theory.  A  modification  in  the  criterion  for 
the  existence  of  a  dip  in  the  output  is  obtained.  It  is  also 
shown  that  effects  due  to  long  range  forces  are  most  noticeable 
at  long  optical  wavelengths  and  when  there  is  a  large  difference 
between  the  lifetimes  of  the  two  la~er  levels. 

The  experimental  system  for  measuring  the  line  profile  of 
a  dc  excited  argon  laser  has  been  designed,  constructed,  and 
made  operative.  Preliminary  data  on  the  "Lamb  Dip"  of  an  argon 
ion  laser  has  been  obtained  as  a  function  of  pressure  and 
excitation.  The  results  of  the  experimental  portion  of  this 
program  and  attempts  to  interpret  the  experiments  in  terms  of 
this  theory  will  be  reported  in  later  reports. 
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INTRODUCTION 

Tho*  objective  of  this  program  is  to  conduct  an 
experimental  and  theoretical  investigation  of  laser  line 
profiles.  The  theoretical  analysis  will  be  based  on  a 
discussion  of  many  particle  distribution  functions  which 
provide  a  single  unified  theory  for  treating  both  the  long 
range  nature  of  the  coulomb  forces  and  the  quantum  mechanical 
interaction  of  the  two-level  system  with  the  optical  field. 
The  effects  of  such  long  range  forces  on  the  steady  state 
conditions  in  a  single  mode  gas  laser  will  be  determined. 

The  experimental  program  will  determine  the  homogeneous  line 
width  in  an  argon  ion  laser  and  relate  experimental  effects 
to  the  prediction  from  the  theoretical  portion  of  the  inves¬ 
tigation. 


KINETIC  THEORY  FOR  AN  ION  LASER 

Introduction 

A  satisfactory  description  of  the  operating  conditions 
for  a  gas  laser  including  effects  due  to  rectilinear  particle 
motion  is  presently  available. ^  A  modification  of  this  theory 
to  include  the  effect  of  short  range  collisions  has  also  been 
developed.'  '  However,  a  theoretical  framework  permitting  the 
systematic  calculation  of  corrections  associated  with  particle 
interactionj  has  yet  to  be  provLded.  The  purpose  of  this  work 
is  to  show  that  such  refinements  appear  quite  naturally  when 
the  dynamics  of  the  laser  medium  is  described  in  terms  of  the 
Wigner  distribution  function. 

Following  the  technique  which  has  been  employed  recently 
in  both  classical  and  quantum  plasmas,  the  Liouville  equation 
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for  an  N-particle  Wigner  distribution  function  is  fiuc  obtained. 
This  distribution  function  contains  coordinates  describing  both 
the  external  center  of  mass  motion  of  each  particle  (and  is 
treated  classically)  as  well  as  the  internal  motion  which  is 
quantized.  From  the  Liouville  equation,  an  equation  for  a  one 
particle  distribution  function  is  obtained  which  contains  an 
integral  over  a  two  particle  function.  This  integral  contains  all 
effects  associated  with  particle  interactions  and  is  exact.  To 
treat  the  long  range  coulomb  forces  that  are  present  in  the  ionized 
gas  laser,  the  collision  integral  is  treated  by  methods  developed 
recently  in  classical  plasma  physics.  Implicit  in  this  treatment 
is  the  standard  assumption  of  plasma  physics  that  one  is  dealing 
with  a  completely  ionized  system  so  that  short  range  collisions 
among  neutral  atoms  may  be  neglected.  Unfortunately,  this 
idealization  is  not  realizable  in  practice.  The  effect  of  long 
range  collisions  on  the  steady  state  conditions  of  an  ion  laser 
operating  in  a  single  mode  is  then  derived  from  a  third  ord  ?r 
perturbation  analysis. 


The  N-Particle  Wigner  Distribution  Function 

We  consider  a  cavity  of  volume  V  containing  N  two-level 

systems  interacting  with  an  electromagnetic  field.  Each  particle 

is  assumed  to  have  a  mass  m  and  the  velocity  distributions  of 

particles  in  each  of  the  two  levels  are  assumed  to  be  Maxwellian 

t  h 

at  the  same  temperature  T.  The  1  such  system  has  internal 
coordinates  p^,  q^ ,  and  external  center  of  mass  coordinates  P^,  Q^. 

The  Hamiltonian  operator  governing  the  dynamics  of  the  N 
two-level  systems  is  assumed  to  be  composed  of  four  types  of  energy 
operators.  They  are  the  kinetic  energy  operator  K.  for  the  center 
of  mass  motion  of  the  i  particle,  the  operator  associated  with 
the  internal  energy  of  an  isolated  two-level  system,  V\  the 
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the  potential  energy  of  a  particle  interacting  with  the  electro¬ 
magnetic  field  and  finally,  V-j  the  two-body  potential  expressing 
the  interaction  between  pairs  of  two-level  systems  in  the  cavity. 

In  actual  laser  systems,  the  scattering  will  be  predominantly 
with  particles  other  than  those  in  either  of  the  two  laser  levels. 
Hence,  to  describe  the  scattering,  it  will  be  neces  iry  to  consider 
a  three-level  system.  Since  optical  transitions  to  this  third 
level  will  be  ignored,  only  a  minor  extension  of  the  two- level 
analysis  will  be  required. 

The  Hamiltonian  is,  then, 


?f  *  Z  (Kj  +  H.  +  VjJ  +  Z  Vij 


i-l 


M 


(1) 


The  two-body  potential  affects  the  internal  dynamics  of 
the  two-level  systems  as  well  as  their  translational  motion.  It 
will  be  assumed  that  these  two  effects  are  additive  so  that  one 
may  decompose  the  potential  into  the  form 


I  V.J  --Z  V„  (IQi-QjD  +  I  Wijtai.iQi -Qjl)  ,  (2) 

i,j  i<i  M 


where  is  the  two-body  potential  influencing  the  c enter  of  mass 

motion  and  W.  .  represents  the  effect  of  the  external  f iel-i  of  the 

th  th 

j  particle  upon  the  internal  dynamics  of  the  i  particle. 

The  first  and  third  types  of  terms  in  the  Hamiltonian  are, 

respectively , 


*  .  a2 

‘  ‘  “  2m  aOi2 

(3) 

Vj  =eE(Qi,t)-qt 

(4) 

The  operator  H^  and  the  second  term  in  Eq.  (2)  are  defined  in 
terms  of  the  Schroedinger  equation 


(Hi+Z  W|j)*i(qi,t)  =  it,  *(qj,  t)  . 
i 


(5) 
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The  spatial  dependence  of  may  be  replaced  by  a 

stochastic  time  dependence.  Furthermore,  since  the  time  required 

- 13 

for  a  collision  to  take  place  (10  sec)  is  much  longer  than  the 
period  of  the  optical  field  (LCT***  sec),  Eq.  (5)  may  be  solved  in 
the  adiabatic  approximation  yielding  wave  functions  of  the  form 

'jrlrt  (Qi  ,  t)  =  Xin^9|)  exp£-iEnt/t>  "  ( •  *  *  A  E  n  (♦')]  ,  n  :  0,b  ,  (6) 


where 


AE„(0 


A En  i Q |  Qn) 


(7) 


and  the  x  and  are  defined  by  the  eigenvalue  equation 

H|Xjn:EnXin.  (8) 


The  consta  ts  a  and  b  refer,  respectively,  to  the  upper  and  lower 
levels  of  an  individual  two-level  system.  It  is  assumed  that  the 
En  contain  an  imaginary  part  which  accounts  for  decay  to  the  ground 
state  labelled  by  the  subscript  £,  i.e. , 

En  =  En-  rh  yn/2  .  (9) 


As  a  result  of  this  damping,  the  number  of  particles  in  states  a  and 
b  decreases  with  increasing  time.  It  will  be  assumed  that  these 
excited  states  are  replenished  by  an  external  pumping  mechanism  to 
be  introduced  later. 

Associated  with  the  total  energy  operator  ^  is  an  N-particle 
wave  function  satisfying  the  Schroedinger  equation 

(*>7,0  =  t  (10) 

where  17  represents  the  set  of  all  position  coordinates  17  =  (Q,q)  *» 

(Q^  > • • *Qn>  qi,...q^).  The  N-particle  Wigner  distribution  function 
may  now  be  introduced  through  the  definition 


Vf.7?.  »)  =  / 


d6Nx 


ix-f 


(27 r)' 


,6N 


* 

V  ( 


.77 


y  X,  1)'P(‘>7«-y  X, 


t) 


(11) 
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where  £  represents  the  set  of  all  momentum  coordinates  £  =  (P,p) 
p-p...p^),  and  the  6N  dimensional  transform  variable  X 
has  components  (T^,...T^,  r^,...  r^)  corresponding  to  the  3N 

external  and  3N  internal  coordinates  respectively. 


The  Polarization  of  the  Medium  and  Reduced  Distribution  Functions 


kf  a  result  of  using  the  Wigner  distribution  function, 
quantum  statistical  averages  are  obtained  from  expressions  closely 
resembling  those  used  in  classical  statistical  mechanics.  In 
particular,  the  polarization  density  of  the  medium  due  to  the  i1"*1 
two-level  system  is 

(?i(Qi,t)  =  / d^O  /d3Np  /d3Nq  /d3Np (-eq,)fn  (e.rj.t)  (12) 


s  ‘  v/d3p|d3p,d3qie|,f,^i»^,»t^  ’ 


(13) 


where  is  the  reduced  distribution  function  obtained  by  integrating 

1  t"h 

fN  over  all  phase  space  except  that  of  the  i  particle,  i.e.  , 


Mti.iji.o  *  v/d6nd6w-%(N(£,^.) 


(14) 


For  later  use,  a  two-particle  reduced  distribution  function  for  a 
pair  of  particles  i  and  j  :*?  defined  in  a  similar  way  as 

s  V  2y’d6(N‘2)^  d6tN_2)T7  f N(c,->7,  t)  -  (15) 

The  one-particle  Wigner  distribution  function  also  has  a 
representation  in  terms  of  one -particle  wave  functions  as 


fi  U\,V\^  =  e_i*  \  +  xi.t)  t  (16) 
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where 

+  l(7?l»t)  =  I  C*) ^(0| ,  •)  Xn  <Q|  ) .  (17) 

The  sum  on  n  extends  over  a,  b  and  c  while  that  on  v is  actually 
over  the  continuous  values  of  translational  momentum.  With  the 

expansion  given  in  Eq.  (17),  the  one-particle  Wigner  distribution 
function  may  be  written 


=  Ftnm<  pi  tQ t » t )  f|mn( p, ,  q ,  ,t )  ,  m,n  =  ^  (18) 


where 


(20) 


From  these  definitions  it  follows  that 


Fjnm  =  F|mn 

f  -  f  * 

'Inm  ”  Mmn 


(21) 


Now  (P  ^(Q^jt)  is  tne  polarization  density  due  to  the  ith 
particle  at  point  Qt  in  the  cavity.  The  polarization  density  at 

some  point  Q  in  the  cavity  due  to  the  entire  assemblage  of  N  two- 
level  systems  is 


<?{  Q,t)  = 


N 

Z<P. 

1=1 


vQi,t)|  =  -n(0)e/d3Pd3pd3qqf,(P,Q,p,q#t)  ,  (22) 

Ri  =0 


where  n^  ^  is  the  sum  of  the  average  number  density  of  particles 

in  all  levels,  i.e.,  n^°^=  n  n,  n 

a  b  '  c  • 

With  the  form  of  f^  given  in  Eq,  (18)  ,  the  macroscopic 
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polarization  density  becomes 

(P(Q, ♦)  5  -  n^j8b0  ^d3PF,0b (P,Q,t)  +  S*c<3 / d3PFlob(P,Q,t)+  8°Cb / d3PFl0b(PlQ,r)j+ c.c.,  (23) 
where 

=#/ d^^IntqJqXntq)  =  ^nm  •  (2<+) 

In  obtaining  Eq.  (23)  if  has  been  assumed  that  the  off  diagonal 
matrix  elements  are  real  and  hence  equal.  Also,  it  is  assumed  that 
there  are  no  permanent  dipole  moments  so  that  the  diagonal  matrix 
elements  vanish. 

The  optical  field  and  the  associated  polarization  of  the 
medium  may  be  expanded  in  terms  of  the  eigenmodes  of  the  cavity 
Un(Q)  which  satisfy  a  vector  Helmholtz  equation. One  then  has 

E<Q,0=lEl>(t)m,(0> 

n  (z->) 

n 

Assuming  that  the  optical  field  in  the  cavity  varies  in  only 
one  direction  (labeled  by  a  Z  axis)  and  that  the  length  of  the 
cavity  in  this  direction  is  L,  a  set  of  eigenfunctions  adequate  for 
the  present  work  are 


U„(Q)=e  sinknZ  ,  (26) 

where  e  is  a  unit  vector  perpendicular  to  the  Z  axis, 

kn=n7r/l_,  (27) 

and  n  is  a  large  positive  integer.  It  'ill  also  be  assumed  that 
the  polarization  of  the  individual  two-level  system  takes  place 
along  e  so  the  q  =  e|q|  and  SP  mn  -  e|c^mn|*  Finally,  since  we  are 
only  concerned  with  the  contribution  of  level  £  to  the  particle 
dynamics  and  not  to  laser  action  between  levels  a  and  c  or  b  and  c, 
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the  dipole  moments  q3  ac  and  ^ will  be  set  equal  to  zero. 

From  Eqs.  (23)  and  (25)  we  then  have 

(P„(t)  =  -8*n(0)~-^LdZUft(Z)/ d5PFl0fr(P,Q,t  )+c.c.  ,  (28) 

where  “  Yab  “  ^ba- 

As  has  been  shown, the  steady  state  operating  conditions 
for  the  laser  may  be  obtained  once  (Pn(t)  is  known  as  a  function  of 
the  cavity  field  beyond  the  linear  approximation. 


Quantum  Mechanical  Liouville  Equation 

The  time  derivative  of  Eq.  (11)  along  with  the  Schroedinger 
equation  (10)  lead  to 


dU 

at 


i  6  N 


(  2rr) 


■m  * 


[jV(i7+-|-X)-^(i7--|-X]V*(T7-Y-X,t)^(77+“X,t)  .  (29) 


The  contribution  to  atN/at  which  is  made  by  each  of  the  four  types 
of  terms  in  ^  will  now  be  considered  separately. 

The  kinetic  energy  term  given  in  Eq.  (13)  is  readily  shown ^ 
to  contribute  the  convective  term 


dh 


dt 


I  f-» 

=  -— -Ip, 

m  ' 


dU 


ao. 


(30) 


The  contribution  from  the  internal  dynamics  will  merely  be 
rewritten  in  the  form 


atN  1  r  d6NX 
at  n"  'h  J  (2 irf 


/ttV d6VK€0X[l  ( Hj  +Z  W(j)4  fMiCl^:t),  (31) 


where  the  subscripts  ±  refer  to  the  displaced  arguments  of  the  type 
introduced  in  Eq,  (29)  ,  and  use  has  been  made  of  the  inverse  Fourier 
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transform  relation  associated  with  Fq.  (11). 
of  fN  resulting  from  the  third  term  in  ^  is 


me  rcice  or 


vuaugc 


dt 


6N 

=  ie£  E(Qi»»)-/^^eiXCTl^i7~Xft)^(^+-L \ft) 


=  e  £  E  ( Q , ,? )  • 

dpj 


(32) 


d  f 


N 


at 


Finally,  the  part  of  the  Hamiltonian  containing  the  inter 
action  due  to  translational  motion  leads  to 

6N 

®‘iir/(27 r)6N  f6  l(f5)  Xi^[UiJ(,Qi+'|-Ti-Qj--~Tj|) 

(33) 

-U„(IQ,-|-T|-Qj+r.Tj|)j  fH(^Vt\)  . 

The  center  of  mass  motion  will  be  treated  classically.  This  limit 
is  reaaily  obtained  by  expanding  the  integrand  in  Fa  (33)  • 

“•••*"* ~  -u ».  t:  i:,zz:zz‘ 

on.  ...over.  ,1, 


dfN 

.  y  /  3  U,j 

d 

3Uj, 

d  x 

dt 

dQ, 

dPj 

dOj 

'  dPj  / 

(34) 


tJ  T  ;  y'  "N  3180  ch*nses  due  to  an  external  numping 

ZTlnT  rePleniSheS  the  excited  abates  a  and  b.  This  change 

term**  X  „  n  "T  "  *  3  Phonological  source 

,.'V'  .  ln  he  quantum  hiouville  equation  satisfied  by  f  , 
The  combination  of  Eqs.  (30).  (31)  (32)  anH  ru  \  a  u  W 

Fields  the  quantum  LiouvUle’equ^^li  “e  ££  ““  ^ 

1  A  (35) 

-I(H,*Zvy1)).]fN({lVi-eIE(Q„i,).£lN  .j/aua  J_.  £Uj)  J_\f  s 

“»  aPn  dQj  dP,  dQ,  dpj  . 
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Kinetic  Equation  Satisfied  by  the  One-Particle 
Wigner  Distribution  Function 

It  has  been  shown  previously  that  the  physically  interesting 
quantity,  the  macroscopic  polarization  density,  may  be  expressed 
in  terms  of  a  one-particle  Wigner  distribution  function.  One 
obtains  the  equation  determining  this  function  by  integrating 
Eq.  ( J5)  over  the  phase  space  of  all  particles  except  one. 

Employing  the  definitions  given  *.n  Eqs.  (19)  and  (20)  this  procedure 
yields 


dft  Pi  df.  .  d  f. 

"77~  +  in  +£  (' wnm+>n(n)  ^  ‘  T~ 

ot  m  dQ|  n,m*o,b,c  dP 


df 


(36) 


where 


M 


6(N-I)  6(N-I)  & 


and 


Xmn=Xnm  =  (  Y  m  +  yn  )  /  2 


(37) 


(38) 


It  will  prove  useful  to  rewrite  the  term  containing  the  external 
field  by  using  the  identity 


d6X 


dp. 


(2tt) 
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where  d6X,  =  d3T,  d3r,  .  In  terms  of  the  functions  F-^nm  introduced 
previously , 


eE(0„t)  ^-=-ieE(0„l)/^e‘iT'  FIM,x"(q,- ^t,)  X„(q,+ \ r,i  •  <40> 


Modification  of  the  Collision  Term  for  Long  Range  Forces 

The  one  term  in  Eq.  (36)  requiring  further  elucidation  is 
the  collision  integral.  Since  the  distribution  function  f ^ 
satisfies  an  equation  containing  integrals  over  f^  which  describes 
three  particle  effects  etc.  ,  a  means  of  closing  this  hierarchy  must 
be  introduced.  Long-range  forces  may  be  treated  by  truncating 
Eq.  (36)  in  a  manner  similar  to  that  employed  in  plasma  kinetic 
theory.  Neglecting  particle  correlation  entirely  and  writing 

f2(£l,£2,'>7l,‘»72,t)  s  fi  (Ci  »^|i » t)  f  |  (f2,T?2,t)  ,  (41) 

a  pair  of  Vlasov-like  equations  may  be  derived  for  the  distribution 
functions  Faa  and  F^  introduced  previously.  It  should  be  emphasized 
that  no  expansion  in  any  plasma  parameter  is  employed;  spatial 
inhomogeneity  in  the  assemblage  of  two-level  systems  is  assumed  to 
be  due  solely  to  the  presence  of  the  externally  applied  optical 
field.  One  finds  that  the  introduction  of  Eq.  (40)  into  Eq.  (36) 
yields 


*1.  Ji 

dt  m 


aoi 


+  X  ( i  ^nm-^  /nm) 


F  f 

1  nm  1  mn 


eE 


|lL+9£.|lL=5,(P|iQ|,t),(42) 


where 


»  €  ■  -"(0)  ~  jV p2d3  02  U|2  [  F00  +  Fbb  +  Fcc  ] 


Introducing  a  potential  through  the  definition 

ff  -  _i± 

3Q| 


:^3> 


(44) 
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and  employing  the  orthogonality  properties  of  the  f. . ,  namely. 

(27rT,)3/d3pd3qf0(p,q)fob(p,q)  =Sjo8(i.  _ 

one  finds 


(45) 


jfqq.  +  Pi  dFaa 
Of  m  0q( 

,  P|  OFqj, 
dt  m  ’  IqT"  +  1 


OFc, 

at" 


+  y.F.  o-i 

n 

'  (Fab-Fbo) 

d<t> 
-  e— - 

do. 

aFao 

Oh  =  Xo 

(46) 

i^abPab 

j  E  •  a  b 

( Foa  -  Fbb)  - 

d<f> 

eac^ ' 

aFab  _ 
ap, 

(47) 

a  fcc 

dQ,  ' 6 

a  <p 

a  Fee 

m 

dQi 

ap, 

-  =  0 

(48) 

where 


Xo(P1,Ql,t)=(2rrh)3/dap,d3qif,coa,(5,7)„l)  . 


09) 


When  these  equations  are  expanded  to  third  order  in  the  optical 
^engt  ,  it  is  found  that  the  equations  governing  F  (°) 
same  as  those  in  the  collisionless  theory. 


and  Fab<'1'*  are  the 
namely , 


dFal0)  P  dFe{0) 


dt  + 

rn  0Q  +  7a  Foa  s  Xo(P,Q,f) 

(50) 

dfob1"  P 

dt  m 

A F  <•> 
0hob 

*  ao 

(51) 

dPoo2>  p 

aF  ^ 

wraa 

dt  +  nn 

ao 

(52) 
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^ob3)  t  P 
dt  +  m 


aF0b(3) 

dQ 


+  iiiobFoi3) 


(53) 


where 


One  also  obtains 


&ob  “  W0b  “  1  Xob  • 


(54) 


IEhl.’  +  £l . +  ya  F0ia-  ,  .  lEs!  =  JL  hi.  f  FobC)  .  Fboll)  ]  (55) 

dt  m  dQi  dQi  d  P  i  ti  L  •* 


13)  .  (3) 

dFob  P|  oFab 

+ - : - +  ifta 


at 


m  dQ 


_  (3)  iE-  frab  f_  (2)  _  (2)  ]  a*  aFob 
bFob  =  - - -  F00  -  Fbt  +  e  — - — -  ,  (->o) 

*  1  J  ao,  ap, 


(2) 

It  is  seen  that  Faa '  '  is  goverened  by  an  inhomogeneous  Vlasov 
equation  and  hence  may  be  obtained  through  the  use  of  standard 
integral  transform  techniques. 


A  Review  of  Collisionless  Results 

Before  considering  the  effect  of  collisions,  the  collisionless 
results  of  Ref.  (1)  will  first  be  recovered.  From  Eq.  (50)  one 
obtains 


F.rv.o.o = / 

*  o 


00 


dre  Yo  X0  (P,Q- Pr/m,t- r) 


(57) 


Assuming  tu"t 

(0) 

x0(P,Q,t)  =  A0$  (P)  ,  (58) 

where  A  a  is  constant  or  varies  slowly  compared  to  the  other  length 
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and  time  scales  of  interest  in  the  problem,  one  may  write 


(00) 


F.r-(P,Q,t)  =  Ao4>t0)(P)/X0  . 

For  use  in  Eq.  (51),  it  is  useful  to  write 


(59) 


.  (00)  (00)  (0)  ,  (0) 
'oa  “^bb  =<**  N/n 


where 


(60) 


N 


Aa  A  b 

77  '  77 


(61) 


and  is  referred  to<i>  as  the  "excitation  density".  The  equation 
obtained  by  interchanging  a  and  b  in  Eq.  (51) 


_  is 

FabO,(P,O,t)-Fba(0)(PlQlt)  =  i8(N/n0))<l>(0>(P) /°dr  (e"ia°bT+e“il2baT) 

•O 

x  U„(0-Pr/m)An(f-r)  , 

where  the  electromagnetic  field  has  been  written 


(62) 


E  (Q,t)  s  eEnAn(t)u„(Q) 


and 


(63) 


=  E„  8VTi  . 


(64) 


Obtaining  a  similar  solution  for  Faa^20^  -  F^20)  from  Eq.  (52) 
and  substituting  the  result  into  Eq.  (53)  leads  to 


Fob  (P,Q,t)  =  (18;  (N/n  )  $l0)(P)  Z  /“dTdrdrV^7"^' 

c=a,b 

X  (e  '^ab7  +e  ,ftba  T  )  An{  t  -  t)  An(  t  -  T  -  r‘)  An(t  -  r  -  T*  -  T") 

X  Un(Q-  Pr/m)Un[o  -  (Pr  +  Pr') /m]  U„[q-  (Pr  Pr'  -8-  Pr") /m] . 


(65) 
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From  Eq.  (28),  the  first-order  contribution  to  the 
polarization  is 


/P Co)  (0)  2  C  3  ft  do) 

(To  (t)  =-n  ^ J  d  Pjo  dZUn(Z)Fob  (P,Q,t)  +  c.c  .  (66) 

Introducing  Fab(10),  carrying  out  the  spatial  integration  and 
neglecting  the  second  harmonic  contribution  yields 


(Pn'0)(t)  =-18  8'N/d3P<J>(0)(P)/Odre‘!n°bTAn(t-r)  cos  k„Pr/m  + 


c.c. 


(67) 


=  - i 8  S* N  A  (i-t) 

"  o 


+  c.c. 


where 


ft0  5  knp0/2m  . 


(68) 


Setting 


A„(l)  =  cos  ( wnt  +  4>) 


(69) 


and  using  the  rotating  wave  approximation^-^  one  obtains 


(P n'0)(t)  =  -  «',(Wn,4^z  (-^U  c.c. 

2  kn  i/0  \  kn  i/0  / 


(/0) 


where 


Z  (C)  *  i  /o°°dte^,-,*/4  =  Zr  +  iZj 


(71) 
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In  Ref.  (1),  it  is  shown  that  the  steady  state  operating 
conditions  fellow  from  a  knowledge  of  that  part  of  (Pn  that  varies 
as  sin  ( t  +<£  ) , 

Writing 

(Pn(t)  =  Cncos(a>nt+ *)  +  Snsin(a.nt +  $)  (72) 

there  results 

2 

S  (,0)  =  -  g  NEn  2j  (  U)n~nob  \  (73) 

"  fiknv0  '  knv0  '  ’ 


The  third-order  contribution  is  obtained  from  Eqs.  (28)  and 
(65)  by  a  similar  but  somewhat  more  lengthy  procedure.  The  spatial 
integration  is  first  carried  out  by  employing  the  formula 


4 


~  /LdZUn(2)U„(Z-a,)Un(Z-a2)U„(Z-  a3) 

L 

[cos  (at  +  a2-a3)  +  cos  (a2 +  a3-<X|)  +  cos  (a3 +  ai- a2)]  . 


(74) 


In  the  present  case 


at  +  a2  -  a3  =  knP(r-r")/m 

a2+a3-ai  =  k„  P(r+2r‘ +  r")/m  (75) 

a3 +ai  -  a2  =  knP  (r +r")/m 

As  has  been  pointed  out,^-^  the  second  and  third  of  these 
terms  give  contributions  with  higher  inverse  powers  of  knvQ  and 
hence  in  the  "extreme  Doppler  limit"  may  be  neglected.  Carrying 
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out  the  momentum  integration  then  yields 


(Pnl30l<t)=4S'(i8»5NE  /“< iT  dr'dr“e'ifl»“r->''V'  n'(T  r"11 
*  c=oVo 


(76) 


x(e‘'^obT '  +e'i^b°r  )An(t-r)  An(t-r-r')  An(t-  r-r'-r")  +c.c.  . 


When  the  three  time  factors  are  written  in  exponential  form,  it 
is  found  that  only  the  terms  (l/6)e'i( +  f  *  +  ‘wn(r+r")  antj 

(1/8)  e~  1  ("n» ♦vj  +  ;  give  rige  t0  denominators  that  will 

have  a  resonance  when  the  various  time  integrations  in  Eq.  (76) 
are  performed.  Retaining  only  these  terms,  Eq.  (76)  becomes 


<P 


(30) 


5Ne-U„„t.f) 


dr  dr'dr"  e'^° 


-ftjV-r")2 


x  e 


•  '  wwr 


1 


(77) 


-g  V  ('  8)5Ne-i<‘V*'*l(_L+  J_  )  [x  UiolW,  ,u,)+  I(«0,w„-w*  )]+c.c.. 


where 


CD,  =  cu 
I  n 


■a 


ob 


( 78) 


and 


Kfl01W,w')sfdT()TVn5(r'',l!+i‘'lHW,r' 

*  r\ 


(79) 
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One  integration  in  this  expression  is  readily  performed  after 
introducing  the  transformation 


x=t,+t2  (80) 

y=ri-r2. 

The  remaining  integrals  ere  expressible  in  terms  of  the  function 
introduced  in  Eq.  (71).  One  finds 

i(n^^')-[z(-^-)+2(1^-)]/[2flo(uW)].  <81> 

Since  Z(-z*)  =  -Z*(z) ,  it  follows  that 

IUl0,aj|lwl)  =  (knVoWl)',Z(cu,/2n0),  (82) 

and 

)  =  (knv0Imw,)  'z }(w(/2ft0).  (83) 

In  terms  of  these  results,  the  third-order  polarization  becomes 


"  32knvc 


-i-  )  e1  ( 2  (uj)  /  knVo)  +  ( IfTvaj|f 1 Z ,  ( OJ,  /  knv0)  ]  (84) 


From  which  one  finds 


(30) 

Sn  = 


tf4E  3N 

8h  3/q  /bknVo 


(W|  /knv0)  [l  +  ybb2^C(wn-Woii)] . 


'Oc'i 


where 


I 

w2+Xob2 


(86) 
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Now  the  operating  condition  for  a  iaser  is  obtained  from 
'  steady  state  solution  of  the  equation 


En  -  En*  (3 n F 


where 


(37) 


(JOB i 

'  Qi.  7  2  U  0  /  E 


(i) 


_IL 

En 


(88) 


*•4  (-£■)-¥ 


(3) 


(89) 


of'free^iLe"6  Q  ""  “Vity  “de  and  is  the  P^ittivity 

obtainedhh  threSh°ld  excicaci°n  density  for  laser  operation  is 
obtained  by  setting  a.  -  0  in  Eq.  (8b).  Neglecting  effects  that 

*  *  °( *  \*  °ne  flnds  that:  the  threshold  density  at  resonance 

*  Wob  )  IS 


Nt= 


*o V0 

~Qn 


(90) 


(90)  -°  ^  E<1'  (87)  ^  introducinS  Eqs.  (73),  (85)  and 

(90),  the  equilibrium  field  amplitude  is  given  by 


2E  2  (wn-«ab)2/(fenVo)2 


1 


yay b  l+Xab^T(u;n-w0b) 


(91) 


wnere 


V  =  N/Nt 


(92) 
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In  obtaining  this  result  the  approximation 

Zj(x)=7rl/2  e"  *2  ,  (93y 

has  bet'.n  used  in  the  numerator  to  yield  the  Doppler  profile.  In 
the  denominator  this  exponential  has  been  replaced  by  unity  since 
the  dip  in  the  response  curve  occurs  in  the  region  iwn-wob|<yob 
for  which  the  exponential  is  essentially  constant  in  the  extreme 
Doppler  limit. 


Solution  of  Vlasov  Equation  for  F 

n  aa 

The  solution  of  Eq.  (55)  may  be  carried  out  by  introducing 
the  transforms 

Poo121  (P.0,0  =  I  e'VQ  /Hj  e"fo1,2l(P,s)  ,  (94) 

*12l(Q,t)  *  •*'  ♦fill  ,  (95) 

h»,2’(P,0,')  «  ^‘‘kiafSr  •*,H(1<|)(P,s)  ,  (96) 


where  h0  ( Pt  Q  ,  t )  is  the  inhomogeneous  term  in  Eq .  (55). 
The  transformed  version  of  Eq.  (55)  is 


(s+r 


a 


„  .  .  vf  (2)  . 

P|/m)  Ja£  +  i  ek£ 


dF, 


(0) 


00 


dP, 


9jt 


(2) 


=  H 


(2) 


(97) 
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From  £q.  (43)  it  follows  that 


(98) 


Equation  (97),  along  with  a  similar  equation  for  FK/(2)plus  a 

homogeneous  equation  for  F  (since  so  *  so  -  m  / 

.  ,  ,  .  cl  vaAUt'e  ^ac  •  hr  “  °)  may  now 

b'  solved  simultaneously.  One  finds 


foTlP.s)  .  #J"(Pi, ,  +  WV  ng«(.) 

"(r;  Us) 


i  k 


k  2 


_a»(0) 

dP 


s  +  Xa  +  ik^P/m 


,  (99) 


where 


and 


r(2°)/D  .  Ha  (p,s) 

Pa,  (P,S)  =  - - -  , 

s  +  Xa  +  i  k^  P/m 


(100) 


«rW"V 


(101) 


Also,  «^(s)  Is  defined  by 


f^(s)  =  €0j[(s)  +  tb£(s)  +  €C/(s)  -  2 


where 


(102) 


eng($)  =  I ~ 


4^/d3p 


ik. 


d<t> 


(0) 


ap 


s  +y'n>  i^’ P/m  *  n  "  °>b,c  .  (103) 
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Third-Order  Polarization  for  the  Ion  Laser 

The  third-order  polarization  follows  in  the  usual  way  from 
Eq.  (28), 


(P  n3)(t)  =  -n^^ob  —  j[Lo  2o't1(Z)  /  d3PF0P(P,Q,t)  +  c.c.  .  (104) 

From  Eq.  (56)  the  appropriate  form  for  Fab^  is  found  to  be 
Fab(3)( P, 0,0  s  ( iE0' ?ob/'h)  O "  r)  U„ (Q - P r/m )  [ Foo2) ( P,  Q  -  Pr/m ,  t  -  r) 

(105) 

-  Fbb2)(P,  Q-Pr/m,t-r)]  +  dr#-|fW£  (2)  (Q-  Pr/m  ,t-r)  ~  F0bl,)  (P,Q-Pr/m,t-r) 

The  spatial  integration  required  in  the  evaluation  of  the  polariza¬ 
tion  density  is 


|-/oLdZUn(Z)U„(Z-VT)e^z  =  y  (2»iocos  knVr-a(?i2nei'iVT-8/,-2ne‘i,,-AVT)  .  (106) 

It  is  found  that  the  dominant  correction  to  the  polarization  comes 

from  the  X  »  0  term  in  Faa^^-F^^\  All  other  terms  yield 

contributions  that  are  higher  order  in  (k  v  )”^.  Hence  we  set 

'no' 


A^(3)=  (Pnl3)-(Pn30)  =  -  fobUEo  SWF>) /odre  in<sbTA„(t-r)  J d3P / 

x  -lk.Pr/m  mn*  / _ j _ _ J _ \ 

kz2(s2+a,p2)  \»  +  y.  +  »kx-P/m  "  s  +  yb  +  Jk^P/m  / 


dS  »(t*r) 

1WT9 

cos  kn  Pr/m 


(107) 


2  (O)  2 

where  ujp  =  4  n  n  e  /m . 

The  is  »  0  form  of  the  dielectric  constant,  namely 


<0($)  =  l  +  s2/wp2 


(108; 


-23- 


F-920479-2 


has  already  been  introduced  in  Eq.  (107)  since  it  is  insensitive 
to  the  limiting  process  required  in  evaluating  this  expression 
at  k ji  -  0. 

Performing  the  s-inversion  first,  the  correction  to  the 
third-order  polarization  may  be  written 

£(p(3)  2  w0  +  Wt+c.c.  ,  (109) 


where 


Wa 


=  /d“T,-|n"“r  A„(t-r) / d!  C06  k„  PT/ 

\  1i  /  2  i  k. 2  Jo  J 


m 


i  k. 


d<S> 

dP 


(0) 


(110) 


iwp+  ik/P/m 


+  (similar  term  with  wp  ■—  -cup)  , 
and  ujp2c  -  4  v  n^'e2/m . 

The  seconc  order  density  variation  to  be  used  in  Eq.  (110) 
is  readily  obtained  from  Eqs.  (50)  -  (52).  The  result  is 


N«°l(t)  =  -^N(iEoT0b/l.)2/odrdT1e"a*2T'2(e'r°T-e_rbT) 
X  (e-i^<JbT'+  e-'^boT')  An  (t  -  r)  An(t  -  r  -  t') 


(HD 


Integration  over  momentum  in  Eq.  (110)  requires  the  evaluation  of 
the  integral 


ik< 


d$ 


(0) 


dP 


e-ik£  P12/mcos  .  pf  /m 


Ya  +  iu>p+  ik^  P/m 


(112) 
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Expanding  the  c'enominator  and  only  retaining  terms  0(k£  *’)  yields 

Ia(t|,t2)=  -  k/e‘(fl°,l)2[l-2(ri0t,)2][l+t2(Xa4-ia;p)]/ (  Xa  +  iwp)  .  (113) 

When  Eqs.  (110)  and  (111)  are  combined  and  the  three  time  factors 
An(t)  are  written  in  exponential  form,  it  is  again  found  that  only 
two  terms  give  rise  to  resonant  contributions  in  subsequent  time 
integrations.  Hence  Eq.  (110)  may  be  written 

W0  =  ^NwPwPa2rob(E0-  S’ob/t*)  T*  ita°b_u,n)T  [l  -  2Ui0r-)2] 

X  (Xa  +  iWp)'2^  dr'  {e,wPT' [l  +  r(/a  +  ia<p)]-e'raT  [l  +  (r+r,)(Xa  +  iwp)]  }  (114) 

x  f  0dr"dT“,en°2 T,"*(e‘r°T“-e-ybT“  )[e*l(flob‘Wn)r"+  e',tabaf  Wn)T'"  1 


+  (similar  term  with  u>p  —  ~up  ) 

After  performing  the  simple  integrations  of  r'  and  r"  ,  the  r,n 
integration  is  given  in  terms  of  the  z  function  introduced  in 
Eq.  (71)  and  the  r  integration  requires 


dCe'S*+i<C(l-2C2)  =  -  y€+0(€2) 

/o°°dCe'C  +i,C  cO-2C2)=-~+0(6) 


(115) 


One  finally  obtains 


VV- 


S3 ob  /  ^-o  “cMjfc  W  I  I 


4/a'  n  Vq)' 


.  ai.> 


The  resulting  change  in  the  polarization  density  may  now  be  obtained 
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from  Eq.  (109'  The  associated  Jn  quadrature  component  is  found 
to  be 


AS. 


(3) 


Fob* 


4yarb*n  v0 


(I^£mZ(i£U)[ 


(117) 


where  L 


kT/4rr  n 


(0)e2 


Finally,  the  equilibrium  field  amplitude  is  given  by 


2  2 

SYbEn 

8t»xaxb 


( wn-ftab)*  /  (kn*'#)* 

I  -  M  e 


I  +  Xa b<5C ( Wn  -  Wob)  +  2(X0-  Xb}2  [(X0knLa)'2  +  (XbknLbr2  ] 


.  (H8) 


Comparison  of  Eq.  (118)  with  the  corresponding  result  of 
Ref. (l)shows  that  the  effect  of  long  range  forces  is  contained 
in  the  last  terra  in  the  denominator  of  Eq.  (118).  This  term  is 
enhanced  by  increasing  the  wavelength  of  the  laser  radiation  and 
by  decreasing  the  Debye  length  associates  with  the  particles  in 
states  a  and  b.  Since  these  Debye  lengths  turn  out  to  be  on  the 
order  of  10”^  cm  while  the  optical  wavelengths  are  10”"*  -  10”^  cm, 
it  is  seen  that  the  correction  term  would  be  quite  small  were  it 
not  for  the  fact  that  the  transition  rate  yb  of  the  lower  level 
may  be  as  much  as  10  yQ  . 

As  a  numerical  example,  consider  l/u.  radiation  in  a  plasma 
having  a  temperature  of  2000°K  and  n(0)  ~  n(K01~5  x  10^  cm”^.  Then 

22  -72  ob 

La  =  =  2  x  10”  cm  .  Assuming  that  y ^  =  10  ya  ,  the  correction 

term  is  0.2.  The  response  of  the  laser  at  frequencies  above  o»ab 
both  with  and  without  this  correction  is  shown  in  Fig.  1.  The 
complete  response  curve  is,  of  course,  symmetric  about  the  frequency 
cuab  .  Although  the  shape  of  the  response  curve  remains  unchanged 
when  long  range  forces  are  included,  a  superposition  of  the  two 
curves  in  Fig.  1  shows  that  the  effect  of  such  forces  is  to  fill 
in  the  center  of  the  curve.  Also,  the  requirement  of  positive 
curvature  at  the  center  of  the  response  curve  yields  as  a  criterion 
for  the  existence  of  a  dip, 
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where 


-|-  =  (ro-Xb)2  [ ( yakni-0 )  2  +  (ybknLb)  2]  _  (120) 

Finally,  Eq.  (118)  shows  that  the  correction  term  vanishes 

if  y  -  y ,  .  This  is  clear  on  physical  grounds  since  in  the 

a  D  <"? )  (2) 
laser  theory  one  finds  that  Faav“'  =  '  when  Ya  and  y^  are 

equal.  These  two  terms  then  have  identical  spatial  dependences 

that  are  180°  out  of  phase  with  each  other.  Since  the  net  force 

acting  on  a  given  two-level  system  is  due  to  the  inhomogeneous 

distribution  of  other  charged  particles,  it  is  clear  that  for 

y  a  ~  the  distribution  of  particles  on  state  a  plus  those  in 

state  b  gives  rise  to  a  completely  homogeneous  spatial  distribution 

of  such  particles  and  no  net  force  is  felt  by  an  individual  two- 

level  system. 
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MEAS'Uk-  .4  NTS  OF  ARGON  ION  LASER  LINE  PROFILE 

Introduction 

In  considering  the  output  power  of  laser  oscillators,  it 
has  been  shown  that  the  saturation  characteristics  of  the  optical 
medium  depend  upon  the  broadened  line  width  of  the  medium  and 
the  frequency  separation  between  modes. v  '  A  representation  of 
this  mechanism  is  most  easily  described  as  a  hole  being  "burned" 
in  the  Doppler-broadened  gain  curve  of  a  gas  laser.  The  "hole 
burning"  refers  to  the  saturation  of  the  atoms  in  the  vicinity  of 
the  frequency  of  the  saturating  radiation. At  lowT  power  levels 
the  distortions  in  the  "hole  width"  are  due  to  the  broadening 
caused  by  collisio  in  a  non-ionized  plasma^)  and  by  the  long 
range  forces  in  an  ionized  plasma (see  the  earlier  section  of 
this  report). 


Experimental  Arrangements  and  Preliminary  Results 

For  a  cavity  of  very  short  length,  the  spacing  of  the 
resonances  may  exceed  the  Doppler-broadened  line  width  and  the 
laser  will  oscillate  in  a  single  mode.  One  expects  oscillation 
at  a  single  frequency  to  cause  gain  saturation  of  the  medium  only 
in  the  vicinity  of  this  frequency.  Two  holes  are  "burned"  in 
the  gain  curve  due  to  the  standing  wave  character  of  the  optical 
fields  in  the  cavity.  If  the  length  of  the  cavity  is  tuned 
through  the  center  of  the  Doppler  profile,  a  "dip"  is  observed 
which  is  due  to  the  fact  that  few7er  atoms  contributing  to  the 
power  output  of  the  laser.  For  the  specific  case  of  the  ionized 
argon  gas  laser,  the  gain  profile  is  so  wide  that  a  short  single 
mode  cavity  could  not  be  conveniently  constructed  and  therefore 
the  "Lamb  Dip"  could  not  be  displayed  by  tuning  the  c/2L  frequency 
through  line  center.  If  the  output  intensity  of  the  argon  laser 
was  extremely  stable  (which  it  is  not),  it  could  be  operated  very 
close  to  threshold  in  one  mode  and  the  "Lamb  Dip"  could  possibly 
be  exhibited;  however,  no  range  of  the  excitation  parameter  would 
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exist . 

An  alternate  solution  which  has  been  successfully  applied 
during  the  reporting  period  is  to  use  a  three  mirror,  variable 
reflectivity  system  at  one  end  of  the  resonator . A  /2  motion 
of  this  system  will  allow  the  two  cavity  system  to  pass  through 
alignment  and  will  provide  a  resonance  for  one  of  the  modes  of 
the  original  cavity.  By  linearly  driving  one  of  the  mirrors  of 
the  variable  reflectivity  cavity  the  resonances  may  be  scanned 
across  the  whole  Doppler-broadened  line  and  will  exhibit  the 
"Lamb  Dip"  as  the  particular  oscillating  modes  are  scanned  in  ~ime. 
The  "Lamb  Dip"  may  now  be  plotted  as  a  function  of  excitation  and 
the  area  of  the  "hole"  represents  the  relation  power  to  be 
achieved  by  a  single  mode  since  the  "hole"  area  represents  the 
homogeneously  broadened  line  width. 

The  experimental  arrangement  is  shown  pictorial] y  in  Fig.  2 
and  schematically  in  Fig.  3a.  A  d.c.  Ar+  discharge  is  contained 
in  a  quartz  capillary  which  is  60  cm  in  length  and  has  a  2  mm  ID. 

The  excitation  can  be  varied  from  3  to  10  amperes.  An  intracavity 
prism  is  incorporated  into  the  design  to  enable  the  various  atomic 
lines  to  be  selected.  The  three  mirror,  variable  reflectivity 
system  is  located  at  one  end  of  the  cavity  and  has  a  multiwafer 
piezoelectric  drive  associated  with  mirror  M,  which  provides  the 
scan  of  the  short  cavity  relative  to  the  long  cavity. 

The  system  has  been  designed,  constructed  and  mad j  operative. 
Preliminary  data,  as  shown  in  Fig.  3b,  has  illustrated  the  "Lamb 
Dip",  and  demonstrated  the  suitability  of  the  technique  for  a 
detailed  investigation  of  the  shape  of  the  hole  burned  in  the 
Doppler  gain  curve.  Detailed  results  of  the  experimental  portions 
of  this  program  and  attempts  to  interpret  the  experim  nts  in  terms 
of  this  theory  will  be  reported  in  a  later  report. 
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SIX-MONTH  STATUS  EVALUATION 

Some  aspects  of  the  theory  of  a  gas  laser  recently  developed 
by  W.  E,  Lamb,  Jr.  have  been  recast  in  a  form  which  more  fully 
displays  the  role  played  by  the  particle  dynamics.  The  Wigner 
distribution  function  has  been  used  to  derive  kinetic  equations 
which  govern  the  external  center  of  mass  motion  of  the  two-level 
systems  as  well  as  their  internal  dynamics.  The  effect  of  long 
range  forces  has  been  discussed  by  treating  the  collision  integral 
in  a  manner  similar  to  that  employed  in  plasma  kinetic  theory.  A 
modification  in  the  criterion  for  the  existence  of  a  dip  in  the 
output  has  been  obtained.  It  has  also  been  shown  that  effects  due 
to  long  range  forces  are  most  noticeable  at  long  optical  wavelengths 
and  when  there  is  a  large  difference  between  the  lifetimes  of  the 
two  laser  levels. 

The  experimental  system  for  measuring  the  line  profile  of 
a  d.c.  excited  argon  laser  has  been  designed,  constructed,  and 
made  operative.  Preliminary  data  on  the  "Lamb  Dip"  has  been 
obtained  as  a  function  of  pressure  and  excitation.  Detailed  results 
of  the  experimental  portions  of  this  program  and  attempts  to  inter¬ 
pret  the  experiments  in  terms  of  this  theory  will  be  reported  in 
a  later  report. 
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